s17aff
s17aff
© Numerical Algorithms Group, 2002.
Purpose
S17AFF Bessel function J (x)
1
Synopsis
[fun_value,ifail] = s17aff(x<,ifail>)
Description
This routine evaluates an approximation to the Bessel Function of
the first kind J (x).
1
Note: J (-x)=-J (x), so the approximation need only consider
1 1
x>=0.
For very large x, it becomes impossible to provide results with
any reasonable accuracy, hence the routine fails. Such arguments
contain insufficient information to determine the phase of
________
/ 2
oscillation of J (x); only the amplitude, / -------, can be
1 \/ (pi)|x|
determined and this is returned on soft failure. The range for
which this occurs is roughly related to the machine precision;
the routine will fail if |x|>~1/machine precision
Parameters
s17aff
Required Input Arguments:
x real
Optional Input Arguments: <Default>
ifail integer -1
Output Arguments:
ifail integer